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Abstract 

The leading-order process for the production of a Z boson and a heavy-quark jet 
at hadron colliders is gQ — > ZQ {Q = c,b). We calculate this cross section at next-to- 
leading order at the Tevatron and the LHC, and compare it with other sources of ZQ 
events. This process is a background to new physics, and can be used to measure the 
heavy-quark distribution function. 
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1 Introduction 

Many signals for new physics at hadron colliders involve electroweak gauge bosons {W^, Z, 7) 
and jets containing heavy quarks (c, b). The prime example is + 4 jets, with one or more 
jets containing a b tag, which led to the discovery of the top quark [HIS El- Other examples 
include signals for the Higgs boson and the superpartners of the known particles IHEj. It 
is therefore crucial to understand the standard-model background from the production of 
electroweak bosons and heavy-quark jets with good accuracy. 

The simplest processes of this type are the production of a single electroweak gauge boson 
and one heavy-quark jet. In the case of the W boson, the leading-order process is gs W~c 
[Hj, which has been calculated at next-to- leading order (with rric nonzero) |7j. For a photon, 
the leading-order process is gQ —>■ •yQ {Q = c,b), which has also been calculated at next-to- 
leading order (with rriQ = 0) |Hlinj- In this paper we consider the analogous leading-order 
process for the Z boson, gQ — >■ ZQ, shown in Fig. ^ which we calculate at next-to-leading 
order (with mg = 0). 




Figure 1: Associated production of a Z boson and a single high-p^ heavy quark {Q = c,b). 

An alternative calculational scheme is to regard gg ZQQ as the leading-order process 
(with rriQ nonzero), and to allow one heavy quark to be emitted coUinear to the beam, 
yielding a ZQ final state. This approach has two drawbacks. First, the expansion parameter 
of this calculation is as In^Mz/mq) rather than as, so perturbation theory is less convergent. 
Using a heavy-quark distribution function sums these coUinear logarithms to all orders, 
resulting in a perturbative expansion in as and 1/ In^Mz/mg) jTUl HH [12] ■ Second, it 
is much more difficult to obtain the next-to-leading-order correction to gg ZQQ than 
to gQ ZQ, both because there is one more particle in the final state, and because 
the heavy-quark mass must be maintained in the calculation of gg ZQQ to regulate 
the coUinear region (one may set mq = in gQ — > ZQ, provided that the heavy-quark 
transverse momentum is much larger than its mass). At present, gg — > ZQQ is known 
at next-to- leading order only for mg = jlHl CH CHI- In contrast, the next-to- leading- 
order calculation of gQ — ^ ZQ (with mq = 0) that we perform in this paper is sufficiently 
straightforward that the next-to-next-to-leading-order calculation may be available in the 
foreseeable future \W\ . 

The process gQ — *■ ZQ is a background to gb —>■ hb, where the Z boson and the Higgs 
boson decay to the same final state {bb, r^r^, or [T71 UHl CH] . In addition, gQ ZQ 

could potentially be used to measure the Q distribution function. The b distribution function 
is needed for the above-mentioned Higgs-production process as well as for inclusive Higgs 
production, bb ^ h [201 |2Il 122] • It is also needed for the single-top production processes 
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qb q't [121123 and gh — > W~t |211I2S1, and the charged-Higgs production process gh H^t 
j2Hl 123 12H] • The process gQ jQ is much more sensitive to the charm distribution function 
than to that of bottom, due to the greater electric charge of the charm quark [201 IHI] • 

At present, the b distribution function is derived perturbatively from the gluon distri- 
bution function fOl 1201 and there is no direct measurement of it. The c distribution 
function is similarly derived from the gluon distribution function, and it agrees well with 
direct measurements. Thus we expect the same to be true of the b distribution function. 
The uncertainty in the b distribution function derives from the uncertainty in the gluon 
distribution function. 

Another source of ZQ events is qq ZQQ, shown in Fig.|2l where either one Q is missed, 
or the two Q's coalesce into a single jet. We show that this process is much more significant 
at the Tevatron than at the LHC, since light-quark distribution functions are relatively more 
important at large values of Bjorken x. This calculation is carried out at leading order with 
a non-zero Q mass. 



fWWZ 



(a) 



Q 
Q 




Figure 2: Representative Feynman diagrams for qq ZQQ. The Z boson may be radiated 
off (a) the initial-state quarks or (b) the final-state quarks. 

In addition, we also calculate Zj production at next-to-leading order, where j denotes a 
light-quark or gluon jet, as shown in Fig. [31 32] [^SliSll • Using a silicon vertex detector to tag 
heavy quarks, the probability that such a jet fakes a heavy-quark jet is around 1%. Taking 
this probability into account, we show that this source of fake ZQ events is comparable to 
genuine ZQ events at the Tevatron, but is relatively less important at the LHC. 



fVWVZ 



9 ^TTOTTO^ 



fWWZ 



(a) 



(b) 



Figure 3: Representative Feynman diagrams for Zj production via (a) gq — * Zq and (b) 
qq Zg. 

All of the above processes can also lead to final states with two jets, with varying numbers 
of heavy quarks. For completeness, we also calculate these cross sections (at leading order). 
When combined with the next-to- leading-order cross section for ZQ, one obtains the inclusive 
cross section for Z plus at least one heavy-quark jet (at next-to-leading order). 
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The next-to-leading-order calculations in this paper were performed with the Monte- 
Carlo code MCFM The leading-order calculations were performed both with this code 
and with MadEvent jHEj- 

2 gQ ^ ZQ at NLO 

The next-to-leading-order (NLO) calculation of gQ ZQ parallels our NLO calculation of 
gb —>■ hb [in], and we refer the reader to that work for a detailed discussion of the calculational 
scheme. The contributing subprocesses are 

• gQ ~^ ZQ (one loop) 

• qQ ZQq 

• gQ ^ ZQg 
•99^ ZQQ 

The subprocess qq — > ZQQ is considered separately in the following section, and is not 
regarded correction to gQ ZQ. 

We work in the simplified ACOT scheme [TIH ITT| . which allows one to neglect the Q 
mass throughout. This is a good approximation, and simplifies the calculation. The error 
made by this approximation is proportional to l/ln^Mz/mq) x rriQ/p^ jTH]. We use the 
dip ole- subtract ion method as formulated in Ref. jSH] to isolate and subtract coUinear 
divergences. 

In our NLO calculation of gQ — > ZQ, we demand one and only one jet with transverse 
momentum pt > 15 GeV within a rapidity range \ri\ < 2 at the Fermilab Tevatron. This jet 
must contain a heavy quark. At the CERN Large Hadron Collider (LHC), the heavy-quark 
jet must have pt > 15 GeV and \ri\ < 2.5. 

If two partons lie within a cone of radius AR < 0.7, we merge them into a single jet 
with four-momentum equal to the sum of the two partons' four-momenta. This is done 
before the cuts described above are apphed to the jets. The NLO process gg — > ZQQ yields 
two heavy quarks in the final state. If they are merged into a single jet, we record it as a 
double-heavy-quark jet. This is only about 1% of the ZQ cross section. 

Similarly, the NLO process QQ' QQ'Z, and processes related by crossing, yield two 
heavy quarks in the final state. However, these processes amount to a correction of less than 
1%, so we neglect them ^19j . 

We list in the first column of Tables ^ and |21 the LO (in parentheses) and NLO cross 
sections for gb —>■ Zb and gc Zc at the Tevatron and the LHC. The Zc cross section exceeds 
that of Zb by 70% at the Tevatron and 35% at the LHC because the charm distribution 
function is larger than that of bottom. This is partially compensated by the fact that the Z 
has stronger coupling to bottom than to charm. The ratio of the Zb and Zc partonic cross 
sections is proportional to 
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More importantly, a silicon vertex detector (SVX) can tag a b jet with an efficiency around 
60%, and a c jet with an efficiency of about 15%. Thus the majority of ZQ events tagged 
with an SVX come from Zb. 

The NLO processes that contribute to ZQ also give rise to final states with more than 
one jet, and we list these cross sections in the remaining columns. These events are grouped 
in three classes. The second column, labeled Z{QQ), corresponds to events with a single jet 
that contains two heavy quarks. As mentioned earlier, this is only about 1% of the ZQ cross 
section. The third column, labeled ZQj, corresponds to events with two jets, one of which 
contains a heavy quark. This is about 1/5 of the ZQ cross section at the Tevatron, and 
about 1/2 at the LHC (for pt > 15 GeV). The fourth column, labeled ZQQ, corresponds 
to events with two jets, both of which contain heavy quarks. This is significantly less than 
ZQj. The ZQQ events arise from the NLO process gg ZQQ (as do the Z{QQ) events). 
The final column sums these various processes, to give the inclusive cross section for Z plus 
at least one heavy-quark jet at next-to-leading order. 

We show in Figs. EHZI the scale dependence of the inclusive cross section for Zb and 
Zc production at the Tevatron and the LHC at leading order (LO) and next-to-leading 
order (NLO). The dependence on the renormalization and factorization scales are shown 
separately. Both scale dependencies are reduced at NLO in comparison with LO. We use 
fiR = Mz, fJ'F = Mz as our default value of the renormalization and factorization scales. 

We also estimate the uncertainty in the NLO inclusive cross section. The ffist uncertainty 
is due to varying the renormalization scale between half and twice its default value of /ir = 
Mz- The second uncertainty, obtained in the same manner, is due to the factorization 
scale. The third uncertainty is from the parton distribution functions [22] • There is also an 
additional uncertainty of 4% due to Sas{Mz) = 0.002. 

We show in Figs. IHHHl the transverse-momentum distributions of both the Z and the 
highest-pT b jet at the Tevatron and the LHC . At leading order these distributions are 
identical, since the Z recoils against the b jet. At next-to- leading order the distributions for 
the inclusive cross section are slightly different, since the Z can recoil against two jets. The 
distributions for Zc production are qualitatively similar. 

3 qq^ ZQQ at LO 

Another contribution to ZQ production comes from qq ZQQ {q = u,d,s), shown in 
Fig. 121 where one Q is outside the coverage of the detector. The dominant contribution 
to the cross section comes from diagrams in which the Z is radiated from the initial-state 
quarks while the heavy quarks arise from gluon splitting, as shown in Fig.|21^a). We maintain 
the heavy-quark mass throughout the calculation in order to regulate the singularity that 
would arise from a gluon splitting to massless coUinear quarks. Although the NLO cross 
section for this process with massless quarks is available [121 CH HH] , the cross section with 
massive quarks is known only at LO. 

The LO cross section for ZQ production via qq —>■ ZQQ is given in the ffist column of 
Tables H and 121 At the Tevatron, the cross section is about 1/3 of gQ — > ZQ. At the LHC, 
the contribution from qq —>■ ZQQ is relatively much less, only about 1/30 of gQ —>■ ZQ. 
This reflects the fact that this process is initiated by a qq collision. At the Tevatron, where 
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Table 1: Cross sections (pb) for Z-boson production in association with heavy-quark jets at 
the Tevatron {^/s = 1.96 TeV pp). A jet hes in the range pt > 15 GeV and |?7| < 2. Two 
final-state partons are merged into a single jet if Rjj < 0.7. No branching ratios or tagging 
efficiencies are included. The labels on the columns have the following meaning: ZQ = 
exactly one jet, which contains a heavy quark; Z{QQ) = exactly one jet, which contains two 
heavy quarks; ZQj = exactly two jets, one of which contains a heavy quark; ZQQ = exactly 
two jets, both of which contain a heavy quark. For the last set of processes, the labels mean: 
Zj = exactly one jet, which does not contain a heavy quark; Zjj = exactly two jets, neither 
of which contain a heavy quark. For ZQ and Zj, both the leading-order (in parentheses) and 
next-to-leading-order cross sections are given. The last column is the next-to-leading-order 
inclusive cross section, which is the sum of four previous columns. The CTEQ6M parton 
distribution functions are used throughout, except for the LO cross sections in parentheses, 
where CTEQ6L1 is used [22]. The factorization and renormalization scales are chosen as 
/^F = /^R = Mz- The uncertainties are from the variation of the renormalization scale, the 
factorization scale, and the parton distribution functions, respectively. 



Cross sections (pb) 


Tevatron 




ZQ Z{QQ) ZQj ZQQ 


ZQ inclusive 


gb^ Zb 
qq Zbb 


(8.23) 10.4 0.169 2.19 0.631 
3.32 1.92 - 1.59 


13.4 ±0.9 ±0.8 ±0.8 
6.83 


gc — > Zc 
qq Zee 


(11.3) 16.5 0.130 3.22 0.49 
5.66 6.45 - 1.70 


20.3t||±0.ll|:^ 
13.8 




Zj Zj] 


Zj inclusive 


m Zg,gq Zq 


(876) 870 137 


ir|in+44 +9 +7 
iUiU _4Q _2 _12 



the typical values of Bjorken x are relatively large, the valence quark distribution functions 
are significant. In contrast, the typical values of x are relatively small at the LHC, so the 
valence quark distribution functions are less important. 

The second column in Tables ^ and |21 lists the contribution from qq ZQQ when the Q 
and Q merge into a single jet. This is comparable in size to the contribution of this process 
to ZQ production. This is due to the aforementioned enhancement that arises when the 
heavy quarks are collinear. This process is the dominant source of events where the Z is 
accompanied by a single jet that contains two heavy quarks. 

The last column gives the contribution of qq —>■ ZQQ when both jets are within the 
coverage of the detector. The Zbb cross section is about 1/2 of the Zb cross section from this 
process at both machines, and the Zee cross section is about 1/3 of the Ze cross section. 
This is in contrast with the gQ — *■ ZQ process, where the ZQQ final state (that arises at 
NLO) is much less than ZQ. 
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Table 2: Same as Tabled except at the LHC {^/s = 14 TeV pp). A jet lies in the range 
Pt > 15 GeV and |?7| < 2.5. 



Cross sections (pb) 


LHC 




ZQ Z{QQ) ZQj ZQQ 


ZQ inclusive 


gb Zb 
qq Zbb 


(826) 649 11.3 304 78.1 
24.3 13.5 - 11.4 


in/in+70+70 +30 

iU^U _go _ioo -50 

49.2 


gc — > Zc 
qq — i> Zee 


(989) 921 8.8 396 61.5 
36.7 41.7 - 11.3 


1390 ± lOOlfo -m 
89.7 




Zj Zjj 


Zj inclusive 


QQ Zg,gQ Zq 


(13500) 11600 4270 


T c;S7n +900 +60 +300 
iOO/U _goo -300 -500 



4 qq ^ Zg, gq Zq at NLO 

A light-quark or gluon jet can fake a heavy-quark jet. With a silicon vertex detector, the 
mistagging rate is typically around 1%. Since the cross section for Zj production is much 
greater than that for ZQ, it is potentially a large source of fake ZQ events. 

We list in Tables Hand 121 the LO (in parentheses) and NLO cross sections for Zj produc- 
tion. The LO processes are qq ^ Zg {q = u, d, s, c, b) and gq ^ Zq {q = u, d, s), as shown 
in Fig. El The NLO processes also contribute to the Zjj final state, and we list that cross 
section in the Tables as well. The Zj cross section is almost two orders of magnitude greater 
than that of gQ —>■ ZQ at the Tevatron. Taking into account the 1% mistagging rate, we 
see that Zj is a significant source of fake ZQ events at the Tevatron. In contrast, at the 
LHC the Zj cross section is not nearly as significant. Thus there will be relatively fewer 
mistagged events at the LHC. Qualitatively similar results are obtained for the ratio of Zjj 
to ZQj. 

5 Conclusions 

The dominant contribution to Zb production is gb —>■ Zb. At the Tevatron, qq — > Zbb also 
makes a significant contribution. Combining these two contributions, the total inclusive 
cross section for Zb production at the Tevatron is about 20 pb. Thus there are about 2000 
inclusive Zb events produced within the coverage of the detector at the Tevatron for every 
100 pb~^ of integrated luminosity delivered. The cleanest decay mode of the Z boson is to 
leptons {i = e, /i), with a branching ratio of 6.7%. Including a 6-tagging efficiency of 60% 
yields about 80 tagged {Z £^£^)b events for every 100 pb^^. There are slightly more 
tagged events when one accounts for the fact that some events contain two heavy quarks, 
either within the same jet or in separate jets. With data sets between 4000 and 8000 pb~^ 
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expected in Run II at the Tevatron, there will be between 3200 and 6400 tagged {Z — > £^£^)b 
events. 

At the LHC, the contribution from qq —>■ Zbb is much less significant, so most of the cross 
section comes from gb — >■ Zb. The total inclusive cross section for Zb production is about 
1090 pb, a factor of about 50 larger than at the Tevatron. Thus there will be an enormous 
number of Zb events at the LHC. 

The total inclusive cross section for Zc production at the Tevatron is about 70% greater 
than that of Zb, while at the LHC, Zc is about 35% greater than Zb. However, the SVX 
tagging efficiency for charm is about 1/4 of the 6-tagging efficiency, so there will be fewer 
tagged Zc events than tagged Zb events. In contrast, the number of 7c events far outweighs 
the number of 76 events at both machines, since charm has twice the electric charge of 
bottom and the charm distribution function is larger than that of bottom. Even taking the 
greater tagging efficiency of bottom into account, the number of tagged 7c events is larger 
than the number of tagged 76 events [HTH IHT|. 

Zj events, where the jet is mistagged as a heavy quark, are a significant source of fake 
ZQ events at the Tevatron, but much less so at the LHC. Thus a larger fraction of the tagged 
ZQ events at the LHC will be from genuine heavy-quark production. For this reason, and 
also due to the relatively smaller contribution from qq — »• ZQQ, the LHC provides a cleaner 
environment for the extraction of the heavy-quark distribution functions via gQ ZQ. 

In addition to the decay Z — £^i~, there will be many Z ^ vv events [BR = 20%), 
which will yield heavy-quark monojets. This will also yield dijet events with large missing 
transverse momentum, with one or both jets containing heavy quarks. 

Regardless of the Z decay mode, the majority of Z + 2 jet events with a single heavy- 
quark tag at the LHC come from ZQj, not ZQQ. At the Tevatron, where qq — > ZQQ is 
relatively more important, the ZQj and ZQQ final states are comparable in size. 
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pp^Zb @ Tevatron 
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Figure 4: Cross section for gb Zb at the Tevatron vs. the renormahzation scale (sohd 
curves) and factorization scale (dashed curves). The ratio of the cross section at scale to 
the cross section at scale fi = Mz is plotted vs. the ratio of the scales. The next-to-leading- 
order (NLO) inclusive cross section (bold) is less sensitive to the scales than the leading-order 
(LO) cross section (regular). 
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pp^Zb @ LHC 
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Figure 5: Same as Fig. HI but at the LHC. 
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Figure 6: Same as Fig. |1[ but for Zc production. 
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Figure 7: Same as Fig. IHl but at the LHC. 
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Figure 8: The transverse-momentum distribution of gb Zb at the Tevatron. Quahtatively 
similar results are obtained for gc — > Zc. The LO curve is the pt of either the Z or the b jet. 
At NLO, the distributions of both the Z and the highest-pT & jet are shown for the inclusive 
cross section. 
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Figure 9: Same as Fig. but at the LHC. 
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